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Abstract 

HLA-A, -B and -DRBl allele frequencies and their haplotype frequencies in 21,918 Chinese residents living in Liaoning 
Province, who were registered as volunteer donors of China Marrow Donor Registry, were investigated. They are composed 
of 93.37% Han Chinese, 5.1% Manchus, 0.57% IVlongols, 0.46% Hui persons, 0.29% Koreans and 0.14% Xibe ethnic group. In 
total eighteen different HLA-A alleles, forty-eight different HLA-B alleles and fourteen different HLA-DRBl alleles have been 
identified. Their frequencies are in agreement with the Hardy-Weinberg equilibrium. For Han Chinese in Liaoning, 1,534 
different HLA-A-B-DRBl haplotypes were identified, with a frequency of higher than 0.01%. A*30-B*13-DRB1*07, A*02-B*46- 
DRB1*09 and A*02-B*13-DRB1*12 are the most frequent haplotypes among Liaoning Han. While Liaoning Han, Liaoning 
Manchu, Liaoning IVlongol, Liaoning Hui and Liaoning Korean share the northern Han characteristic haplotypes, all minority 
ethnic groups with the exception of Liaoning IVlanchu have developed their own unique HLA profiles. This dataset 
characterizes the HLA allele and haplotype frequencies in the Liaoning area and suggests that it is different from those in 
other parts of China and ethnic groups, which implicates transplant donor searching strategies and studies on population 
genetics. 
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Introduction 

Human Leukocyte Antigens (HLA) are located on the 
chromosome 6 in tlie 6P21.31 region and are the most 
polymorphic human genes with more than 8,000 alleles [1]. 
HLA genes play an important role in organ and Hematopoietic 
Stem Cell (HSC) transplantations. In order to provide stem cell 
transplant patients with more chances to find HLA-matched 
unrelated donors, China Marrow Donor Program (CMDP) 
launched the nationwide HLA typing program for marrow donor 
registry volunteers recruited over 23 different provinces in China 
in 2001. Our laboratory,contracted for HLA typing with CMDP, 
has typed about 22,000 volunteers from Liaoning Province for 
HLA-A, -B, -DRBl loci from 2003 to 2006, and we uploaded the 
related data to the CMDP database annually. Up to now, CMDP 
has enlisted 1 .8 million potential HSC donors and facilitated more 
than 3,800 HSC donations, including 127 donors from our 
Liaoning bank. 

In earlier studies, we reported single locus HLA allele 
frequencies with a smaller sample size [2—4]. Here we present 
the largest HLA-A, -B and -DRBl typing dataset in Liaoning 
Province. The samples were collected from 14 different cities 
throughout Liaoning Province. Although the entire samples were 
not randomly selected since they were exclusively potential donors 



for the China marrow donor bank, they contain a lot of useful 
information for analyzing the HLA gene diversity in Liaoning and 
its relationship with other parts of China. 

Northern Han Chinese constitutes the main ethnic group of 
Liaoning, although many minority ethnic groups have lived here 
for thousand years. The Liaoning region was actually the home of 
over five millions of Manchu population who ruled China from 
1616 to 1912 during the Qing Dynasty. There is also significant 
presence of Mongolian, Uyghur, Korean, Hui and Xibe ethnic 
populations in Liaoning. Our aim was to analyze the HLA allele 
and haplotype frequencies among Liaoning residents and compare 
them with other northern Chinese, southern Chinese, and 
minorities living outside Liaoning. The results are reported as 
below. 

Materials and Methods 

Ethics Statement 

This project was approved by Medical Ethics Committee of 
Liaoning Blood Center, China. AH participants provided their 
written informed consent and all clinical investigations were 
conducted according to the principles expressed in the Declaration 
of Helsinki. 
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Characterization of Study Subjects 

The volunteers were age 18-50 year old residents of the 
Liaoning, China who were randomly recruited through donor 
drives organized by local Red Cross societies in college campus, 
local community centers and shopping centers. The volunteers 
were multi-ethnic and consisted of Han Chinese, Manchu, 
Mongolians, Hui, Koreans, and Xibe men and women. The 
male/female ratio was 53.75%/46.25% in the study subjects. 

Samples and DNA Extraction 

For each volunteer a sample of 5 ml whole blood was collected, 

divided equally, and stored at — 80°C until processed. The DNA 
was extracted by using the mini column whole blood DNA 
extraction kit from Tiangen (Beijing, China). The DNA concen- 
tration was 40-80 ng/|il and optical density at 260/280 was 1.6- 
1.8. 

HLA Typing 

All individuals were typed as HLA class I (HLA-A, and -B) and 
HLA class II (HLA-DRBl). For the Sequence Specific Primers 
(SSP) method, exons 1-4 of HLA-A and B loci and exon 2 of 
DRBl locus were covered. For the Sequence Specific Oligonu- 
cleotide Probe (SSO) method, exons 2 & 3 of class I loci and exon 
2 of class II locus were typed. Amplification was accomplished on 
the GeneAmp PCR system 9700 (Applied Biosystems, Foster City, 
CA, USA). 

About a third of the samples were typed with the SSP method 
by means of the commercial kits from Life Technologies 
(Carlsberg, California, USA) and Biotest (Biotest AG, Dreieich, 
Germany). About two thirds of samples were typed with the SSO 
method by use of LABType HLA-A,-B,-DRB1 kits from One 
Lambda (Canoga Park, California, USA). The SSP typing results 
were analyzed with Unimatch plus 3.0 and Biotest HLA-SSP 
typing software provided by the manufactures, while the SSO 
results were analyzed by using HLATools software provided by 
One Lambda. When ambiguous results were produced from the 
SSP and SSO typing methods, the Sequencing Based Typing 
(SBT) kit -SBTexcellerator (Genome Diagnostics B.V., Utrecht, 
Netherlands) was applied to r(;sol\-e unambiguous results. For 
HLA-A and -B, exons 2 to 4 were sequenced, and for DRB 1 only 
exon 2 was sequenced. The HLA-SBT typing results were 
analyzed by means of SBTengine software 3.0 (Genome Diag- 
nostics B.V., Utrecht, Netherlands). For the purpose of revolving 
ambiguous results at the serological level, the high-resolution 
group-specific SSP kits from Olerup SSP (Sweden) were deployed. 
Detail procedures were followed strictiy according to the User 
Instructions provided by the kit manufacturers. 

Statistics Metliods 

Allele frequencies of HLA-A, -B and -DRB 1 were calculated 
with the direct counting method. The deviation from Hardy- 
Weinberg equilibrium was calculated by test [5] . The haplotype 
frequencies were analyzed with the Arlequin software 3., 5 [6] based 
on the Expectation-Maximization (EM) algorithm [7]. The 
Linkage Disequilibrium (LD) of three-locus haplotype was 
calculated using the SPSS 13.0 software package according to 
the report [8], and the LD values, ranged from +1 to -1, reflect 
linkage disequilibrium intensity. The genetic distances between 
different populations were calculated as previously described by 
Nei [9] and a phylogenetic tree was constructed based on the 
allelic frequencies of HLA-A, -B and -DRB 1 and the Neighbor- 
joining method with the PHYLIP-3.695 package. 



Consolidation of HLA Typing Results 

The majority of DNA samples were typed at the intermediate 
resolution by the SSP and SSO methods. In order to consolidate 
the data analysis, all SSP, SSO and SBT data were converted into 
the first field resolution with the exception of B14, B15, B40 
groups which were divided into corresponding serological groups. 

Results 

Ethnic Background of Volunteers 

Liaoning is a northeast province consisting of 84% Han Chinese 
and 16% minority Chinese. The location of Liaoning and its 
geographic neighbors are shown in Figure 1. The volunteers 
consisted of 93.37% Han Chinese, 5.1% Manchus, 0.57% 
Mongolians, 0.46% Hui persons, 0.29% Koreans and 0.14% 
Xibe ethnic persons (see Table 1). Overall Han Chinese is 
overrepresented by about 10%, while the minority groups 
are underrepresented by about 50% relatively to their population 
size. 

Frequencies of HLA-A, -B and - DRBl Alleles in Liaoning 
Han Chinese 

Among the 21,918 donors, 18 alleles for HLA-A, 48 alleles for 
HLA-B and 14 alleles for HLA-DRBl locus were detected. Their 
frequencies are listed in Table 2. The distribution of the allele groups 
are in concordance with the Hardy-Weinberg equilibriitm (p>0.05). 
In Liaoning Han Chinese, the 3 most frequent allele groups for each 
locus are respectively A*02, A*24, A*ll, B*13, B*15(62), B*40(61), 
DRB1*15, DRB1*09 and DRB1*12. B15 is the most diverse HLA- 
antigen family in Liaoning Han Chinese. In total, 9 different 
serological groups within B15 family have been detected, 
including B62(7.239%), B75(4.342%), B71(1.928%), B63(0.303%), 
B72(0.137%), B15(0.464%), B77(0.083%), B76(0.068%) and 
B70(0.032%). The following HLA alleles were found ver\' rarely in 
Liaoning Han Chinese (<0. 1 %): A*25(0.042%), A*34(().() 1 7 "/«), A*66 
(0.090%), A*69(0.078%), A*74(0.015%, B*14(64)(0.08r^), B15*(70) 
(0.032%), B*15(76)(0.068%), B*15(77)(0.083%), B*40(40)(0.015%), 
B*42((0.034%), B*47(0.034%), B*53(0.042%), B*73(0.027%) and 
B*78(0.002%). The foUowing HLA alleles were not detected: A*36, 
A*43, A*80 and B*82. 

Frequencies of HLA-A, -B and - DRBl Alleles of the 
Minority Ethnic Groups 

The Manchu ethnic group constitutes the largest minority 
ethnic group in Liaoning with more than five million residents. 
Their HLA-A, -B, and -DRBl alleles show high similarity to the 
Liaoning Han Chinese. The HLA allele frequencies vary less than 
20% among alleles with a frequency greater than 0.5%. Currently, 
670,000 Mongolians are living in Liaoning and make up 1.53% of 
total Liaoning population. In our study, 126 Mongohan samples 
were included. Their frequencies of A*01, A*68, B*37, B*57 and 
DRB 1*10 are significantiy higher than those in Liaoning Han, 
while those of B*15(62), B*15(75), B*27 and B*48 are much less 
common than those in Liaoning Han. For Hui ethnic group, their 
allele frequencies vary significantiy from those of Liaoning Han. 
A*01, B*07, B*44, B*54, B*57 and B*67 are 50% to 300% less 
common in Liaoning Hui than in Liaoning Han, while their 
B*14(65), B*18, B*38, DRB1*01 are 70% to 250% more common 
in Liaoning Hui than in Liaoning Han. For Koreans living in 
Liaoning, A*03, A*32, B*13 and B*37 have significantiy lower 
frequencies than those in Liaoning Han; on the contrary, the 
opposite is true for B*59, B*67 and DRB 1*04. Only 31 samples 
with Xibe ethnic origin were included in our study subjects. Major 
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Figure 1. The Map of the People's Republic of China. The studied regions are marked with their name box. 
dol:1 0.1 371 /journal.pone.0093082.g001 



variations in HLA allele frequencies between Xibe and the rest 
groups were observed, even so the frequency results of Xibe are 
not reliable due to their small sample size. 



Haplotype Frequencies of Liaoning Han Chinese 

Based on the statistic analysis, 1,534 different HLA-A-B-DRBl 
haplotypes were identified from 20,465 samples with frequencies 
of higher than 0.01%. Twenty-five of diem are classified as 



Table 1. Characterization of 21,918 bone marrow donor volunteers. 



Ethnic Group 


Sample N 


Percentage % 


Liaoning Population 


Percentage % 


Han 


20465 


93.37 


36746911 


84.00 


Manchu 


1118 


5.10 


5385000 


12.31 


Mongol 


126 


0.57 


670000 


1.53 


Hui 


101 


0.46 


264000 


0.60 


Korean 


64 


0.29 


241052 


0.55 


Xibe 


31 


0.14 


133000 


0.30 


Others 


13 


0.06 


306360 


0.70 


Total 


21918 


100 


43746323 


100 


Sex(Male/Female) 


11778/10140 


53.74/46.26 




50.63/49.37 



The ethnic information is based on the registration information filled by the volunteers. The Liaoning population information is based on the 2010 public survey results 
conducted and published by the Liaoning province government (http://www.gov.cn/gzdt/2011-05/13/content_1863235.htm). 
doi:1 0.1 371/journa!.pone.0093082.t001 
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Table 2. Gene frequencies (GF) of HLA-A, -B, and -DRBl in 21,918 bone marrow donor volunteers In Llaoning by etiinic groups. 





Alleles 


Han % (N = 20465) 


Manchu°/o (N = 1118) 


Mongol% (N = 1 26) 


Hul% (N = 101) 


Korean°/o (N = 64) 


Xibe% (N = 31) 


A«01 


4.686 


5.501 


6.746 


1.980 


3.906 


3.226 


A*02 


32.883 


33.274 


31.349 


33.663 


34.375 


27.419 


A«03 


4.630 


3.801 


5.159 


5.941 


1.563 


3.226 


A«11 


15.935 


15.206 


13.095 


10.891 


14.063 


19.355 


A*23 


0.330 


0.134 


0.397 






1.613 


A*24 


16.015 


1 7.039 


15.079 


17.327 


19.531 


8.065 


A»25 


0.042 












A«26 


3.220 


3.309 


4.365 


4.951 


4.688 




A*29 


0.904 


0.581 




0.990 


0.781 


1.613 


A*30 


7.166 


7.021 


5.556 


6.931 


4.688 


14.516 


A^Bl 


4.173 


4.875 


7.143 


5.941 


4.688 


4.839 


A*32 


1.573 


1.521 


1.984 


2.475 


0.781 


4.839 


A«33 


7.286 


6.574 


7.143 


7.921 


10.156 


8.065 


A*34 


0.017 












A«66 


0.090 


0.045 










A*68 


0.955 


0.984 


1.587 


0.990 




3.226 


A«69 


0.078 


0.134 


0.397 




0.781 




A*74 


0.015 












Unknown 


0.002 












B*07 


3.494 


2.996 


2.778 


1.980 


3.906 


3.226 


B*08 


0.909 


0.626 


0.794 


0.495 


0.781 


1.613 


8*13 


12.167 


12.567 


12.698 


12.871 


6.250 


1 7.742 


B»14 


0.002 












6*14(64) 


0.081 


0.045 


0.794 




0.781 




B*14{65) 


0.330 


0.268 




2.475 






8*15 


0.464 


0.537 




0.495 






8*15(62) 


7.239 


7.737 


4.365 


9.406 


8.594 


6.452 


8*15(63) 


0.303 


0.179 


0.794 


0.990 




1.613 


8*15(70) 


0.032 


0.045 










8*15(71) 


1.928 


2.504 


2.381 


1.980 


1.563 


3.226 


8*15(72) 


0.137 


0.134 


0.397 








8*15(75) 


4.342 


4.025 


1.984 


4.951 


3.125 


6.452 


8*15(76) 


0.068 


0.045 










8*15(77) 


0.083 


0.089 










8*18 


0.540 


0.447 


0.794 


2.475 


0.781 




8*27 


1.740 


1.610 


0.397 


0.990 


1.563 


1.613 


8*35 


5.487 


5.456 


8.730 


7.426 


7.813 


4.839 


8*37 


1.610 


1.699 


3.968 


1.485 


0.781 


4.839 


B*38 


2.311 


2.102 


1.984 


6.931 


1.563 




8*39 


1.906 


1.878 


2.778 


2.475 


3.125 




8*40 


0.489 


0.626 


1.190 




0.781 




8*40(40) 


0.015 


0.045 










8*40(60) 


6.355 


7.156 


5.952 


3.960 


4.688 


6.452 


8*40(61) 


7.515 


7.290 


8.730 


4.951 


9.375 


3.226 


8*41 


0.173 


0.045 




0.990 






8*42 


0.034 


0.045 










8*44 


5.089 


4.070 


4.762 


2.475 


4.688 


6.452 


8*45 


0.132 


0.045 










8*46 


6.966 


7.111 


6.746 


6.436 


4.688 


8.065 



8*47 0.034 
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Table 2. Cont. 





Alleles 


Han % (N = 20465) 


Manchu°/o (N = 1118) 


Mongol% (N = 1 26) 


Hul% (N = 101) 


Korean°/o (N = 64) 


Xibe% (N = 31) 


BM8 


3.819 


3.444 




1.984 


2.970 


5.469 




4.839 


BM9 


0.188 


0.089 




0.794 


0.495 








8*50 


0.902 


0.939 




1.190 


1.980 


0.781 




3.226 


B»51 


7.017 


7.111 




4.762 


7.426 


4.688 




6.452 


B»52 


3.655 


3.220 




3.175 


4.951 


3.906 




4.839 


8*53 


0.042 


0.045 














8*54 


3.132 


3.846 




2.778 


0.990 


4.688 






8*55 


1.710 


1.744 




3.571 




2.344 






8*55(54) 


0.024 
















8*56 


0.279 


0.268 














8*57 


1.654 


2.191 




4.365 


0.990 


1.563 






8*58 


4.298 


4.562 




2.778 


3.465 


6.250 




4.839 


8*59 


0.098 


0.089 








3.125 






8*67 


1.033 


1.208 




1.190 


0.495 


2.344 






8*73 


0.027 
















8*78 


0.002 
















8*81 


0.147 


0.045 




0.397 










DR81*01 


2.71 


1.97 




4.37 


6.44 


2.34 




4.84 


DR81*03 


0.03 


0.09 




0.40 










DR81*03{17) 


3.59 


4.03 




2.78 


4.46 


3.13 




4.84 


DR81*04 


10.27 


11.36 




10.71 


15.35 


17.97 




6.45 


DR81*07 


11.66 


12.48 




13.49 


11.88 


7.81 




25.81 


DR81*08 


6.05 


6.57 




4.76 


5.45 


7.81 




3.23 


DR81*09 


13.56 


13.28 




11.90 


10.40 


11.72 




11.29 


DR81*10 


1.60 


1.43 




3.57 


0.99 


0.78 




3.23 


DR81*11 


6.64 


5.50 




6.35 


6.44 


7.03 




8.06 


DR81*12 


12.00 


12.84 




11.90 


12.87 


6.25 




4.84 


DR81*13 


6.31 


5.23 




4.76 


3.96 


7.81 




4.84 


DR81*14 


6.64 


6.31 




7.14 


4.46 


9.38 




12.90 


DR81*15 


17.37 


17.31 




15.87 


15.84 


17.97 




9.68 


DR81*16 


1.56 


1.61 




1.98 


1.49 









doi:l 0.1 371 /journal.pone.0093082.t002 



frequent haplotypes based on a frequency of 0.5%. The most 
frequent haplotypes are A*30-B*13-DRB1*07(4.88%), A*02- 
B*46-DRB1*09 (2.38%) and A*02-B*13-DRB1*12 (1.61%). In 
Table 3, 50 most common HLA-A-B-DRBl haplotypes are listed 
for Llaoning Han and other Llaoning ethnic groups. The 
distribution of Llaoning Han haplotypes based on their frequen- 
cies is shown in Table 4. The 7 most frequent haplotypes make up 
about 14% of total haplotype frequencies. It was found that 229 
haplotypes had frequencies of higher than 0.1% and made up 
about 68% of the whole haplotype frequencies. The rare 
haplotypes with frequencies of lower than 0.01% contribute to 
32% of total haplotype frequencies. 

Linkage disequilibrium (LD) is one of the most important 
characteristics of the HLA system. Forty-nine of the 50 
commonest haplotypes in the Llaoning Han population were in 
positive LD (see Table 5). A*30-B*13-DRB1*07, A*02-B*46- 
DRB1*09, A*02-B*13-DRB1*12, A*33-B*58-DRB1*03(17), 
A*33-B*58-DRB1*13, A*02-B*46-DRB1*08, and A*33-B*44- 
DRB1*13 are common haplotypes with the significant linkage 



disequilibrium (both haplotype frequency and LD are higher than 
1%). 

Haplotypes of Different Ethnic Groups 

Focusing on the ten most frequent haplotypes of Liaoning Han 
(marked in bold in Table 3), the Liaoning Manchus share seven of 
them; the Liaoning Mongols share five of them; the Liaoning Hui 
persons share three of them; the Liaoning Koreans share eight and 
Liaoning Xibe persons share only two of them. Each ethnic group 
shows their unique common haplotypes. For Liaoning Koreans, 
A*02-B*15(62)-DRB1*15 is die most common haplotype, but it is 
not shared by any other ethnic groups within the top seven 
haplotypes. Haplotypes specific for Liaoning Mongols include 
A*24-B*40(61)-DRB1*09, A*02-B*35-DRB1*04, A*24-B*40(61)- 
DRB1*15 and A*24-B*52-DRB1*15. Haplotypes only common 
for Liaoning Hui are A*24-B*15(62)-DRB1*04, A*33-B*14(65)- 
DRB1*01, A*02-B*15(75)-DRB1*04 and A*l 1-B*52-DRB1*15. 
The Xibe ethnic group is exceptional because they have several 
common haplotypes not found in the commonest 50 haplotypes of 
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any other Liaoniiig ethnic group. A few examples are A*02-B*51- 
DRB1*11, A*02 -B*40(60)-DRB1*07 and A*02 B*48- DRB1*14. 

Genetic Distance and Phylogenetic Tree 

The phylogenetic tree is constructed with HTA-A, -B and - 
DRBl allele frequencies presented in Table 2 and the published 
data for the northern Han Chinese from the Shaanxi region [10], 
southern Han Chinese from the Guangdong region [1 1], Koreans 
from Seoul Korea [12], Minnan ethnic group from Taiwan [13], 
Thai people in Thailand [14] and Japanese in Japan [15]. 
Liaoning is a province in the northeastern part of China as shown 
in the map of China (Figure 1). Shaanxi is located in northwestern 
China, which was the center of ancient Chinese history. 
Guangdong Province resides at the very southern end of mainland 
China. Korea is adjacent to Liaoning on its southeast side while 
Mongolia is adjacent to Liaoning on its northwest side. Taiwan, 
Thailand and Japan represent other parts of Southeast Asia. 
Genetic distance here measures genetic divergence between 
populations within a species. A smaller genetic distance indicates 
a closer genetic relationship between two populations, whereas a 
larger genetic distance indicates a further genetic relationship [16]. 
The study on genetic distances between Liaoning Han and the rest 
of study groups revealed that Liaoning Manchus have the shortest 
genetic distance to Liaoning Han Chinese (0.002084), followed by 
Liaoning Mongols (0.015929), Shaanxi Han Chinese (0.020565), 
Liaoning Hui persons (0.027006), Liaoning Koreans (0.036637), 
Guangdong Han Chinese (0.065632), Liaoning Xibe persons 
(0.107017), Seoul Koreans (0.108577), Thai (0.151685), Taiwan 
Minnan (0.228235) and Japanese (0.260845). The phylogenetic 
tree is shown in Figure 2. 

Discussion 

The technology of HLA typing has significantly evolved since 
the PCR method was introduced by MuUis and Faloona [17]. 
Nowadays, there are several kinds of PCR-based HLA typing 
being used: DNA amplification with sequence specific primers 
(SSP) [18], sequence specific oligonucleotide probes (SSO) [19], 
single-stranded conformation polymorphism (SSCP) [20], se- 
quence-based typing (SBT) [21], and DNA chip technology 
[22]. These HLA typing techniques are of much higher accuracy 
as well as reliability than the conventional serologic typing 
method, and also may facilitate the standardization of HLA 
tyjaing processes [23]. The SSP and SSO methods are low- 
intermediate resolution techniques at two-digit level. The advan- 
tage of the SSO method based on Luminex system lies in its high 
throughput, because 96 samples can be completed in one day by a 
skilled technician. SBT is a gold-standard [24], high-resolution 
method at the four-digit level that was used to resolve 
the ambiguous types and to identify novel alleles. In this paper, 
the data of HLA genotyping from 2003 to 2006 were based on the 
SSP and SSO methods. The SBT method was used only to resolve 
the ambiguous types and to identify novel alleles. In 2007, we 
began to use the SBT method in the routine HLA genotyping. The 
high-resolution data are being analyzed and will be reported in 
another paper. 

HLA functions as a useful tool to clarify the difiFerences between 
genetic makeup of Liaoning residents and other parts of China. 
The Han Chinese is divided into two major groups, namely 
northern Han and southern Han [25] . The HLA allele frequencies 
of the Liaoning Han show close similarity to the frequencies of the 
Shaanxi Northern Han, because both have less than 20% 
variation in allele frequencies. Shaanxi played a central role in 
ancient northern Chinese history for thousands of years. This is 
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Table 4. Distribution of HLA-A-B-DRBl haplotypes in Liaoning Han (N = 20,465). 




Frequency Range 


Number of Unique Haplotypes 


% of Total Haplotypes 


>0.01 


7 


13.62 


0.01-0.001 


222 


54.20 


0.001-0.0001 


1305 


29.65 


Total 


1534 


97.47 



doi:1 0.1 371 /journal.pone.0093082.t004 



predictable given that Liaoning is one of the three northeastern 
Chinese provinces that was populated by the geographic 
expansion of the Northern Han Chinese. 

There are striking similarities between Liaoning Han and Liaoning 
Manchu ethnic group in HLA allele and haplotype frequencies. This 
genetic closeness is also supported by the genetic distance calculation. 
This can be explained by the fact that Manchus are descended from 
the Jurchen people with North East Asia origin (http://en.wildpedia. 
org/wild/Manchu_peopleT^Origins_and_early_history). There is 
difiFerence in Manchu HLA allele frequencies between our study 
and previous one [26] . For example, A*33, B*38, B*46, and DRB* 1 4 
have frequencies of 6.57%, 2.10%, 7.11%, 6.31% in the current 
study, but 2.8%, 0.9%, 4.1% and 12.3% in the previous study 
conducted with Manchu residents living in Harbin. Harbin is located 
about 600 km to the north of Liaoning. It is not clear about what 
contributes to this difference beside difiFerent geographic locations. 
Liaoning was the historic home for Manchu. It is possible that more 
intensive gene exchange took place between Liaoning Han and 
Liaoning Manchu in the Liaoning region. The HLA allele 
frequencies of Liaoning Mongol in our dataset are similar to those 
published previously [27,28]. Both B*37 and B*57 distinguish 
Mongol from the northern Han based on much higher frequencies 
in Mongol. 

The HLA haplotype analysis is widely used in human 
population genetics, anthropological studies, as well as the optimal 
marrow donor bank size planning. It carries more specific 
information than allele frequencies. A number of studies have 
addressed the Chinese HLA allele and haplotype distribution from 
different regions including Sichuan, Jiangsu, Fujian, Guangdong, 
Xi'an, Yunan [29-34]. However, there is no report of any large 
sample study about HLA haplotypes in Liaoning. 

Linkage disequilibrium is the common characteristic of HLA 
genetics. This is true for Liaoning Han population according to 
Table 5. The top four most common haplotypes found among 
Liaoning Han in this study include A*30-B*13-DRB1*07, A*02- 
B*46-DRB1*09, A*02-B*13-DRB1*12 and A*33-B*,58-DRB1*03 
(17) which actually consist of allele groups ranking the 5''^ (A*30, 
B*46) and 10* (B*58, DRB 1*03(1 7)) positions in the allele 
frequency ranking (Table 2). B*57 and B*37 rank at the 20* 
and 2r' places in the frequency hst; but their corresponding 
haplotypes, A*01-B*37-DRB1*10 and A*01-B*,57-DRB1*07 
ranked within the top 10 most frequent haplotypes. The LDs 
are similar to the findings in the studies conducted in Shaanxi and 
Inner Mongolia, China [10,28]. 

The fact that HLA-A-B-DRBl haplotype frequencies of 
Liaoning Han show high similarities to those in Shaanxi can be 
expected since both of them belong to the northern China regions 
(the North of Yangtze River). There is below 30% difference in 
haplotype frequencies among 18 out of the top 20 haplotypes 
when Liaoning is compared with Shaanxi (The data were not 
shown). 



Korea is adjacent to Liaoning on its southeast side. Our 
haplotype data do show that the Korean population is much closer 
to the northern Han than to the southern Han. Nine out of top 10 
haplotype rankings in Liaoning Han are also listed in the top 10 
most frequent haplotypes in Liaoning Korean. Korean do have 
own unique haplotypes that show a very different frequency 
pattern compared to Han. The most common haplotype A*02- 
B*15(62)-DRB1*15 in Liaoning Korean is five-time less frequent 
than that in Liaoning Han (Table 3). A*33-B*58-DRB1*13 
(3.91%) and A*33-B*44-DRB1*13 (3.13%), the very common 
haplotypes in Liaoning Korean, reflect also a 3-fold lower 
frequency in Liaoning Han. 

It is well known that A*30-B*13-DRB1*07 is the most common 
haplotype in the northern Han Chinese, while A*02-B*46- 
DRB1*09 is the most common haplotype in the southern Han 
Chinese [29-34]. Indeed, the A*30-B*13-DRB1*07 turned out to 
be the most common haplotype in all ethnic groups in our study 
with the exception of Liaoning Korean. 

A*02-B*35-DRB1*04 and A*24-B*52-DRB1*15, two high- 
frequency haplotypes appear to be unique to the Mongolian 
population, i.e. both of them clearly present a lower frequency in 
Liaoning and Shaanxi Han, but much less in southern Han (The 
data were not shown). 

Our phylogenetic tree study has placed Liaoning Han, Liaoning 
Manchu, Liaoning Mongol, Liaoning Hui and Liaoning Xibe in 
the same close cluster as expected. In addition to the Liaoning 
Han cluster, Liaoning Hui differs substantially from others. A*01 is 
significandy less common in Liaoning Hui than in Liaoning Han 
(1.980% vs. 4.686%), and B*14(65), B*18 and B*38 are 3-7 times 
more common in Liaoning Hui than in Liaoning Han volunteers. 
We have noted that Liaoning Manchu, Liaoning Hui, Liaoning 
Korean are quite difiFerent in their HLA profiles from the same 
ethnic groups reported previously. The differences between 
Liaoning Manchu and Harbin Manchu were already discussed. 
Liaoning Hui also differs in many alleles from Qinghai Hui [26]. 
A*30, A*31, B*14(65), and B*18 show 2.5-12-fold higher allele 
frequencies in Liaoning Hui than in Qinghai Hui, while A*01, 
B*27, and B*40(61) behave in the opposite way. Liaoning Korean 
display a large genetic distance to Liaoning Han, but it clearly 
differs from the Korean residents in Seoul, Korea [11]. It is 
possible that Liaoning Koreans have undergone substantial gene 
exchange with the Liaoning Han population in the Liaoning 
region and have acquired certain HLA characteristics from the 
Liaoning Han while still maintaining Korean unique HLA gene 
composition. The genetic distance between Shaanxi Han and 
Liaoning Han is three times shorter than that between Guangdong 
Han and Liaoning Han. This explains the fact that both Shaanxi 
and Liaoning belong to northern China while Guangdong 
represents a southern China province. Thailand, Taiwan/Minnan 
and Japan lie relative close to each other in the same branch since 
they are all from the Southeast Asia region. Japan population 
shows the most distinct features with the longest distance from the 
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Table 5. Linkage disequilibrium (LD) of HLA-A-B-DRB1 haplotypes in 


Liaoning Han (N = 20,465). 






Haplotype 


HF 


LD 


1 


A»30 B*13 DRB1«07 


4.88% 


4.7783% 


2 


A*02 B*46 DRB1*09 


2.38% 


2.0694% 


3 


A*02 B*13 DRB1*12 


1.61% 


1.1299% 


4 


A*33 B*58 DRB1*03{17) 


1.31% 


1 .2988% 


5 


A*33 B*58 DRBn3 


1.16% 


1.1402% 


6 


A*a2 BM6 DRB1*08 


1.15% 


1.0114% 


7 


A*33 BM4 DRB1*13 


1.14% 


1.1166% 


8 


A*a2 BMa(61) DRB1*09 


0.91% 


0.5749% 


9 


A*01 B*57 DRB1*07 


0.80% 


0.7910% 


10 


A*01 B»37 DRB1*10 


0.79% 


0.7888% 


11 


A*02 B*15{75) DRB1*09 


0.79% 


0.5964% 


12 


A*02 B»15(62) DRB1*15 


0.78% 


0.3665% 


13 


A*n B*15(62) DRB1*04 


0.73% 


0.6115% 


14 


A*ll B*15(75) DRB1*12 


0.68% 


0.5970% 


15 


A*24 BM0(61) DRB1*09 


0.66% 


0.4968% 


16 


A*33 BM4 DRB1*07 


0.65% 


0.6068% 


17 


A*24 B*54 DRB1*04 


0.64% 


0.5885% 


18 


A»11 B*13 DRB1*12 


0.58% 


0.3473% 


19 


A*02 B*40{61) DRBns 


0.55% 


0.1208% 


20 


A*03 B*07 DRB1*15 


0.55% 


0.5219% 


21 


A*ll B*52 DRBl'lS 


0.53% 


0.4288% 


22 


A»02 BMa(60) DRB1*15 


0.52% 


0.1570% 


23 


A*24 B*40{61) DRB1*15 


0.51% 


0.3009% 


24 


A»02 B*4a(61) DRB1*12 


0.51% 


0.2135% 


25 


A*02 B*51 DRB1«09 


0.51% 


0.1971% 


26 


A*02 B*15(75) DRB1*15 


0.49% 


0.2420% 


27 


A*24 B*15(62) DRB1*04 


0.48% 


0.3609% 


28 


A*n B*13 DRB1*15 


0.46% 


0.3069% 


29 


A*32 B»52 DRB1*15 


0.46% 


0.4598% 


30 


A*02 B*15{75) DRB1*12 


0.46% 


0.2887% 


31 


A*02 BM8 DRB1*09 


0.44% 


0.2443% 


32 


A*03 BM4 DRB1*13 


0.43% 


0.4151% 


33 


A*24 B*13 DRB1*12 


0.42% 


0.1862% 


34 


A*24 B»35 DRB1*15 


0.41% 


0.2574% 


35 


A*02 B*35 DRBns 


0.41% 


0.0966% 


36 


A*11 B*40(60) DRB1*12 


0.40% 


0.2785% 


37 


A*24 B*15{62) DRBl'IS 


0.40% 


0.1986% 


38 


A»a2 B*15(62) DRB1*04 


0.40% 


0.1555% 


39 


A*02 B*15{62) DRBn2 


0.39% 


0.1044% 


40 


A»a2 B*48 DRB1*15 


0.39% 


0.1719% 


41 


A*02 B*13 DRB1*07 


0.38% 


-0.0865% 


42 


A*n B*40{60) DRB1*09 


0.38% 


0.2427% 


43 


A*02 B*51 DRB1*14 


0.37% 


0.2168% 


44 


A*11 B»51 DRB1*09 


0.36% 


0.2084% 


45 


A*02 B*15{71) DRB1*04 


0.36% 


0.2949% 


46 


A*24 BM6 DRB1*09 


0.35% 


0.1987% 


47 


A*ll B*40{60) DRB1*15 


0.35% 


0.1741% 


48 


A*24 BMa(61) DRB1*12 


0.34% 


0.1956% 


49 


A*02 B*50 DRB1*07 


0.34% 


0.3054% 


50 


A*02 BM6 DRB1*12 


0.34% 


0.0651% 


cloi:l 0.1 371 /journal.pone.0093082.t005 
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Guangdong Han 



Japan 



■ TaiwanMinnan 



Thailand 



Korean 



Llaoning Korean 
- Llaoning Hul 



Llaoning Manchu 
Llaoning Han 

Llaoning Mongol 

Llaoning Xlbe 



Shaanxl Han 

0.1 

Figure 2. Phylogenetic Tree constructed by the neighbor-joining method display the relationship between Llaoning Han, other 
Llaoning ethnic groups Included In our study, Guangdong Han, Shaanxl Han and populations from Korea, Thailand, Taiwan 
(Minnan) and Japan. 

doi:1 0.1 371 /journal.pone.0093082.g002 



rest of all the study groups. This can be explained with the fact that 
Japan is most isolated nation from their neighbors during the 
thousand years of history. 

Finally it is important to emphasize that our data may not 
reflect the true ethnic makeup of the Liaoning population due to 
the fact that the minority ethnic groups are less represented in the 
marrow donor registry population. The frequency data are also 
based on local bone marrow donor volunteers who could be a 
biased group in terms of representing the entire Liaoning 
population. We have no way to verify the ethnicity of samples 
from minority groups. 

In summary, we have analyzed HLA-A-B-DRBl allele and 
haplotype frequencies of Liaoning residents based on the local 
marrow donor registry volunteers classified by various ethnic 
backgrounds and compared to the populations from other parts of 
China including its geographic neighbors, Koreans and Mongo- 
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